1 



(N 

On 

o 

Q 

00 



> 

m 

(N 
O 
(N 

On 



X 



Cosmic Strings on the Lattice 

A. K. Bukenov a , A. V. Pochinsky, M. I. Polikarpov b , L. Polley c and U.-J. Wiese d 

a University of Alma-Ata, Kazakhstan 

b ITEP Moscow, 117259, Russia 

c Universitat Oldenburg, 2900 Oldenburg, Germany 

d HLRZ Jiilich, 5170 Jiilich, Germany 

We develop a formalism for the quantization of topologically stable excitations in the 4-dimensional abelian 
lattice gauge theory. The excitations are global and local (Abrikosov-Nielsen-Olesen) strings and monopoles. The 
operators of creation and annihilation of string states are constructed; the string Green functions are represented 
as a path integral over random surfaces. Topological excitations play an important role in the early universe. In 
the broken symmetry phase of the U(l) spin model, closed global cosmic strings arise, while in the Higgs phase 
of the noncompact gauge-Higgs model, local cosmic strings are present. The compact gauge-Higgs model also 
involves monopoles. Then the strings can break if their ends are capped by monopoles. The topology of the 
Euclidean string world sheets are studied by numerical simulations. 



1. INTRODUCTION 

During cooling, the early universe has under- 
gone different phase transitions. Depending on 
the symmetry that is spontaneously broken, var- 
ious excitations, for example, domain walls, cos- 
mic strings or monopoles, arise as topologically 
stable objects Jl|,||. The U(l) symmetric scalar 
fields in the 4-dimensional field theories give rise 
to string- like topological excitations. In a bro- 
ken U(l) symmetry phase, the strings are sta- 
ble solutions of the classical equations of motion. 
One distinguishes two types of cosmic strings: 
global and local ones. Global strings arise when 
a global t/(l) symmetry breaks spontaneously, 
whereas local strings are due to the breakdown of 
the C/(l) gauge symmetry. A well known exam- 
ple is the local Abrikosov-Nielsen-Olesen string 
present in the Higgs phase of the noncompact 
abelian gauge-Higgs model with a local U(l) sym- 
metry j3|Q]. An even simpler example is the 4- 
dimensional global U(l) scalar field theory with 
a field <f> = |<t>| exp(zt,c). In two dimensions, ze- 
ros of the complex scalar field, $(2;) = 0, are 
located at isolated points x. Going around x, the 
phase ip may change by 2 7m, where n € Z is 



the topological characteristic of the vorticity of 
the scalar field. The relevant homotopy group is 
Hi [t/(l)] = *Z. In three dimensions, the zeros of 
the scalar field are located on lines (the global cos- 
mic strings), and in four dimensions, on surfaces 
(the world sheets swept out by the strings during 
their time evolution). The strings are topologi- 
cally stable excitations, i.e., they are insensitive 
to small deformations of the field 

The dynamics of cosmic strings and monopoles 
is usually investigated in semiclassical physics. 
In the early universe, however, quantum effects 
played a crucial role. The classical description 
is inadequate when complicated dynamical sit- 
uations arise in which strings or monopoles are 
condensed. The description of phase transitions 
may also require the inclusion of quantum effects. 
Due to the topological nature of the problem, it 
is essential to formulate it in a nonperturbative 
framework, as provided by the lattice regulariza- 
tion. For global cosmic strings, we start with a 
complex scalar field = |$| 2 exp(itp) with a quar- 



tic potential V($) = A(|$|' 



') , and consider 



the limit as A — > 00. Then only the compact 
variable tp € [— ir, it] remains and the theory is re- 
duced to the global U(l) lattice spin model (the 
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four-dimensional XY model). A cosmic string 
then manifests itself by the vorticity of a pla- 
quette and the corresponding dual plaquette be- 
longs to the Euclidean string world sheet. The 
local cosmic strings results from the gauging of 
the U(l) symmetry. On the lattice, this leads to 
the noncompact and the compact abelian gauge- 
Higgs models. The theory with compact gauge 
fields has monopoles as additional topological ex- 
citations. Since the string theory is formulated 
on the lattice, one has complete nonperturbative 
control of the dynamics. In particular, one can 
study the string tension or the question of string 
condensation in the early high temperature phase 
by standard lattice techniques like numerical sim- 
ulations or strong coupling expansions. 

In the present publication, we only give the ex- 
plicit form of the creation operators of the strings, 
the details of the derivation will be published in 
the subsequent paper. The quantization of the 
global cosmic strings is discussed in ref.||. 

2. GLOBAL STRINGS 

Let us consider the 4-dimensional U(l) lattice 
spin model in the Villain formulation |lj . Its par- 
tition function is given by 

Z = E / + "^exp(-|||d^ + 27r/|| 2 ). (1) 

We use the notations of the calculus of differ- 
ential forms on the lattice ||. Dip denotes the 
integral over all site variables, ip; d is the exte- 
rior differential operator, dtp is the link variable 
constructed as usual in terms of the site angles 
p. The scalar product is defined in the standard 
way, e.g., if tp and ip are the site variables, then 
(tp, ip) = J2 S < f 9 ( s )V'( s )i where i s the sum over 
all sites s. The norm is defined as: ||a|| 2 = (a, a); 
therefore ||d</? + 2ttI\\ 2 implies summation over all 
links. Tliifz^^d) denotes the sum over all con- 
figurations of the integers / attached to the links 
C\. It occurs that the partition function ([!]) can 
be represented as follows [||: 

Z K = exp(-2 7 r 2 K (*fc,A" 1 *fc)). (2) 

*fcG^(C 2 ) 
<5*fc=0 



Where *x denotes the object dual to x, the cod- 
ifferential 8 = *d* satisfies the rule for partial 
integration (ip,8tp) = (dtp pip). The summation 
in (||) is carried over the integer variables k, 
which belong to plaquettes Ci- The condition 
8k — means that the summation is performed 
over the closed two-dimensional objects defined 
by k. The surface elements interact with each 
other via long range forces described by the in- 
verse Laplacian. Physically, these forces are due 
to the massless Goldstone bosons in the broken 
phase of the original spin model. The action of 
the random surface model is nonlocal and differs 
from the Nambu action which is proportional to 
the surface area. It can be expected that the ran- 
dom surface model is equivalent to a lattice the- 
ory of closed strings, and that the closed surfaces 
are actually the string world sheets. To justify 
this statement one must construct creation and 
annihilation operators of string states. 

3. STRING CREATION OPERATORS 

The creation of a cosmic string (as a nonlo- 
cal object) requires the use of nonlocal opera- 
tors. Global (local) cosmic strings are surrounded 
by a cloud of Goldstone (gauge) bosons, just as 
charged particles are surrounded by their photon 
cloud. Creation operators for charged particles 
were first constructed by Dirac , whose idea was 
to compensate the gauge variation of a charged 
field, $(x)' = $(x) exp(ia(x)), by a contribution 
of the gauge field representing the photon cloud: 

$ c (x) = $(x) cxp(i J d 3 yB t (x - y)Mv))> ( 3 ) 

where diBi(x) = 8(x), and Ai(x)' — Ai(x) + 
dia(x) is the photon field. The gauge invariant 
operator $ c (x) creates a scalar charged particle at 
the point x, together with the photon cloud sur- 
rounding it. Our construction of string creation 
operators |J is based on the same idea, and it is 
very similar to the construction of soliton creation 
operators suggested by Frohlich and Marchetti 
||. In fact, it is a generalization of their con- 
struction of monopole sectors in 4-dimcnsional 
U(l) lattice gauge theory. We perform the duality 
transformation of the original theory defined by 
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the partition function ([!]), and obtain the hyper- 
gauge theory of integer valued fields on the dual 
lattice. The Wilson loop W(C), constructed from 
the auxiliary gauge fields, is gauge invariant, but 
not hypergauge invariant. We make W(C) hyper- 
gauge invariant by surrounding it by the cloud of 
the hypergauge field. This hypergauge invariant 
operator can be shown to create the closed string 
on the curve C. In terms of the variables k (g), the 
expectation value of the string creation operator 
has the form: 

< U c>=^k x ( 4 ) 
E exp{-2n 2 n{*{k + D c ), A- X *(fc + D c ))}; 

8*k=-S c 

Dc, being the counterpart of the function Bi{x) 
in eq.(||), depends on the curve C: S^*Dc — 6c, 
since the operator of the creation of the string 
should act at a definite time slice, we use the 
three-dimensional operator of the codifferentia- 
tion 5^ ; 5c is the lattice delta function which is 
equal to unity on the links belonging to the curve 
C, and equal to zero on the other links. It is im- 
portant that because of the condition S*k = —5c, 
the summation in (||) is performed over closed 
surfaces, and those bounded by the curve C. This 
is exactly what one would expect intuitively: a 
string world sheet opens up on curve along which 
the string is created. In the general case, the 
curve C may consist of several closed loops: 5c = 
J2i^Ci', then Dc = Y^i^d- Moreover, if we set 
5c j = — 1 on several closed loops, then the string 
is annihilated on these loops. Placing the loops Ci 
and Cj on different time slices, we may construct 
operators corresponding to string scattering and 
decay processes. 

The expectation value of the string creation op- 
erator in terms of the original fields tp is: 

exp(-/3||d ¥ 3 + 2ir5A- 1 (D c - p c ) + 2ttZ|| 2 ). 

Here the integer valued field pc satisfies the equa- 
tion 5^ 3 \*Dc — * Pc) = 0. It can be shown that 



* pc is the analog of the (invisible) Dirac string. 
The Dirac string connected to the monopole is a 
one-dimensional object, while * pc, being defined 
on the plaquettes, is a two-dimensional object. 

4. LOCAL COSMIC STRINGS 

Local cosmic strings are topological excitations 
in the abelian gauge-Higgs model whose partition 
function in the Villain form is given by: 



z = E m E / v * ( 6 ) 

exp(-/3||d0 + 2imf - ~||dp - 6 + 2tt/|| 2 ). 

If the gauge field is noncompact (the integration 
over 6 is from — oo to +oo), then the excitations 
are closed local strings. Below we discuss a more 
interesting case, when the gauge fields are com- 
pact (— 7r < 9 < 7r) and the monopoles are present 
in addition to strings. It occurs that the strings 
may be open, with monopoles at their ends. The 
general Green function for this case consists of 
the creation and annihilation operators of strings 
and monopoles: 

<n^n^iP->= w 

i j k 

\ E / ™ E / ^ 

n6^(C 2 ) w leK{c x ) J n 
exp(-/3\\d9 + 2?m + 2nSA- 1 (B - w)|| 2 ) 
exp{~\\dip + 2ttI-8 + 2tt5A- 1 (D - p)\\ 2 ). 

It follows from the derivation of this formula 
(which is skipped) that 5 6c t + Ylj — 
J2 k 5 Xk = 0; therefore open strings carry 
monopoles and antimonopoles at their ends. In 
the case of creation of monopoles at points Xj and 
antimonopoles at points Xk, the field B satisfies 
the equation: 5^*B = £\ 5 Xj -J2k^ k ; Dirac 
string *uj 6 Z is such that 5^{*B - *u) = 0; D 
and p are defined by the equations: 5^*D = 
5 C + *B- d {3) (D - p) = B - w. It is easy 
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to show that the Green function (0) is invari- 
ant under the deformation of the Dirac string: 
J = uj + d£, p' = p + C 

5. COSMIC STRING DYNAMICS AND 
NUMERICAL SIMULATIONS 

In the high temperature phase of the early 
universe, the U(l) symmetry was unbroken and 
strings were condensed. This was possible be- 
cause the string tension vanished and string cre- 
ation required no energy. As a consequence, one 
expects that the strings formed clusters perco- 
lating through the universe. After the phase 
transition, the U (1) symmetry gets spontaneously 
broken and the string network freezes. Because 
the string tension no longer vanishes, the strings 
become massive. Small strings shrink and de- 
cay, whereas large cosmic strings survive as sta- 
ble massive structures catalyzing the formation of 
galaxies. Such scenario can be studied within the 
formalism discussed above. 

The string condensate can be calculated nu- 
merically by direct measurement of the expecta- 
tion value < Uc >■ A similar measurement of 
the monopole condensate in the compact electro- 
dynamics was performed in rcfs.]9|,|l0|]. Another 
quantity that can be studied numerically is the 
string tension of the cosmic strings. We are per- 
forming these calculations at present. 

The simplest objects related to the cosmic 
strings are the defects in the four-dimensional 
XY model. In the given configuration of the spins 
<p, each plaquette carries vorticity *j — ^([^i — 

<P2hir + Vf2 - Vzh-K + [<£>3 - <Pihn + [<£>4 ~ <Pl]2w), 

where (fk arc the angles on the corners of the pla- 
quette, and [a]2ir is a mod 2ir. If *j ^ 0, then 
we may ascribe vorticity j to the plaquette dual 
to the original one. The surfaces formed by these 
dislocations are closed. The proof is very simple: 
*j = ^d([d(p] 2w ) = ^d(dp + 2np) = dp, where 
the integer p is such that (dip + 2-xp) G (— tt, 7t]. 
Now, using equalities d = *S* and S 2 = 0, we get 
6j = 0. These world sheets are characterized by 
the Euler number: Ne — n p — m + n s , where n p , 
ni and n s are the number of the plaquettes, the 
links and the sites belonging to the world sheet 
respectively. The number of handles is defined by 



.9=1- N E /2. 

We have studied numerically the topology of 
the world sheets in the four-dimensional XY 
model on the lattices 6 4 — 8 4 for the different val- 
ues of k. Usually, there are several disconnected 
closed world sheets in each configuration of the 
spins. It occurs that below the phase transition 
(k < Kc), in the region of the expected conden- 
sation of the strings, there exist one world sheet 
of the size comparable to that of the lattice; this 
world sheet contains a lot of handles. There are 
also satellites: small disconnected objects with a 
simple topology (no handles). After the phase 
transition (k > Kc), almost all world sheets have 
no handles. It turns out that the number of han- 
dles Qi taken into account with the weight propor- 
tional to the area S% of the corresponding world 
sheet i (i.e. < g >= J2i9i^i/J2i >%) i s the order 
parameter of the system. Numerical simulations 
show that < g for k < kc and < g >= 
for k > kc- 

Two of the authors (MIP and LP) would like 
to thank the HLRZ in Julich for hospitality. The 
work of MIP and AVP has been partially sup- 
ported by the grant of the American Physical So- 
ciety. 

REFERENCES 

1. A. Vilenkin, Phys. Rep. 121 (1985) 263. 

2. R. H. Brandenberger, Int. J. Mod. Phys. A 2 
(1987) 77; J. Phys. G 15 (1989) 1. 

3. A. A. Abrikosov, Zh. Eksp. Teor. Fiz. 32 
(1957) 1442; Sov. Phys. JETP 5 (1957) 1174. 

4. H. B. Nielsen and P. Olesen, Nucl. Phys. B61 
(1973) 45. 

5. M.I. Polikarpov and U.-J. Wiese, Preprint 
HLRZ 90-78, Julich, Germany. 

6. P.Becher and H. Joos, Z. Phys., C15 (1982) 
343. 

7. P.A.M. Dirac, Canad. J. Phys. 33 (1955) 650. 

8. J. Frohlich and P. A. Marchetti, Euro. Phys. 
Lett. 2 (1986) 933; Commun. Math. Phys. 112 
(1987) 343. 

9. L. Polley and U.-J. Wiese, Nucl. Phys. B356 
(1991) 629. 

10. L. Polley, M.I. Polikarpov and U.-J. Wiese, 
Phys. Lett. 253B (1991) 212. 



